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■ Abstract 

A calculation of thermal gluon decay shows that this process contributes significantly 
to strangeness production in a quark-gluon plasma. Our analysis does not support recent 
claims that this is the dominant process. In our calculations we take into account the 
resummed form of the transverse and longitudinal parts of the gluon propagator following 
the Braaten-Pisarski method. Our results are subject to the uncertainty concerning the 
estimate of the damping rate entering the effective gluon propagator. 



1 Introduction 



A possible signal for quark-gluon plasma formation in heavy-ion collisions is the enhancement 
of the production of strange particles. The original proposal by Rafelski and Miiller m was 
followed by extensive discussion in the literature H, ||, H, ^, ^ 0] • In this context thermal gluon 
decay has been recently discussed]^, 01- It has been claimed recently that the process g ^ qq 
dominates for a wide range of quark masses |P, 0. Normally, the gluon cannot decay into a 
strange quark-antiquark pair because its thermal mass is too low. Even for the optimistic case 
where on takes the coupling constant g = 2 in a plasma with two massless quarks, the gluon 
mass is given to lowest order in perturbation theory by 



mg = ^gT. (1.1) 



For a temperature of T= 200 MeV this gives = 267 MeV which is below the threshold for 
the production of strange quarks. The important observation by Altherr and Seibert is that 
in addition to acquiring a thermal mass, gluons also acquire a width, determined by the large 
damping rate, of the order gf^Tp]. This is the reason why thermal gluon decay into a heavy 
quark-antiquark pair is allowed, even though the gluon mass is below the threshold for strange 
pair production. 

In this paper we present a systematic re-evaluation of the production rate of massive quarks 
in a quark-gluon plasma due to the processes of quark-antiquark annihilation, gluon fusion and 
thermal gluon decay in the spirit of Altherr and Seibert. Since the production rate depends 
strongly on the damping rate we take a more conservative approach in its estimation. Our main 
point is that even with the parameters chosen in Ref. p, 0] we cannot support the claim that 
gluon decay is the dominant mechanism for strange quark production in a quark-gluon plasma. 
To the best of our knowledge, the gluon fusion mechanism, originally proposed by Rafelski and 
Miiller [m, remains the leading process. 

The paper is organized as follows. In section 2 we briefly review the properties of thermal- 
gluon propagators and the damping rate. In section 3 we calculate the production rates. In 
section 4 we present results and concluding remarks. 

2 Gluon Propagator and the Damping rate 

The effective gluon propagator at finite temperature in the Feynman gauge is given by 

zDj (go, q) = -tS'^'lP^Mqo, q) + P^^^Liq,. q)i (2.1) 
where Pj^ and P^^, are transverse and longitudinal projectors respectively, and 

AT,L(go, q) = ^2 ( (2.2) 

- UT,L{qo,q) 

where Q"^ = q^ — q^- The real transverse and longitudinal parts of the gluon self energy in the 
high temperature limit are respectively given by 



RenT(go,g) = 




41^1n^ 
q 2g qo-q 



(2.3) 
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and 



ReUL{qo,q) = -ml ( 1 - 



Q Qo-q 



(2.4) 



The positions of the poles in the propagator ( p.2| ) are determined by the dispersion relations 

ql = q' + IlT,LiQo,q). (2.5) 

If a pole is located at 

qo = ^T,L + hT,L (2.6) 

then the imaginary shift of the pole 'Jt,l is related to the imaginary part of the self energy 
through 

7r,L = Res(AT,L)ImnT,L, (2.7) 
where Res(A) is the residue of the propagator given by 



Res(A^ 



T,L 



or, explicitly in terms of ujt^l 



Res(Ar) 



-1 



Thus we can write 
A(Q) 



as 



Res(Ai)-i 



Q2 - Ren 



dqo 



-UJT+ + 2 2 : 

-UJL + — + ^. 



(2.8) 



+ 



iRes(A)-i7 



(Q2 _ Ren)2 + Res(A)-272 (Q2 _ YieUy + Res(A)-272 ' 



(2.9) 
(2.10) 

(2.11) 



where we suppressed subscripts T, L, i. e. A, 11 and 7 are either transverse or longitudinal. 
This expression will be used to replace the mass-shell 5-function for thermal gluons. 

The imaginary part of the pole in (|2.6|) gives the damping rate of the plasma oscillations 
1^, |TD| In the following we consider transverse gluons and we set 7 = 7^. The damping rate is 
related to the so called gluon magnetic mass rrimag, or the inverse magnetic screening length at 
high temperature. Unfortunately, the exact relation between 7 and rrimag is not known. A closed 
expression for the damping rate has been derived by Pisarski in the limit where rrimag ^7 0: 
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where the thermal gluon mass is given by 



(iVc 



2 



g V- . ■ 2 ' 9 

The magnetic mass at high temperature is of the form 



mag 



cn9 T 



(2.12) 



(2.13) 



(2.14) 
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where is a number depending on the gauge group and cannot be calculated by a perturbation 
expansion. Lattice estimates |Tl|, |T^ for SU(2) 



C2 = 0.27 ±0.03 



(2.15) 



have been confirmed by recent semiclassical calculations [|T3|. So far, no reliable estimate exists 
for SU(3) H^. The best one can do is extrapolate the SU(2) value by [pTSj] 



C3 



Expanding the log in ( |2.12|) in powers of 
finds 

l = {l + riV^ 



2 



C2 



(2.16) 
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mag and retaining only the leading terms one 
2 \ 



In 



1.09681... 



mag , 



(2.17) 



where = if we keep the leading log term only. If the next-to-leading term is included then 

N 



V 



Aire 



N 



(2.18) 



To check the consistency of Pisarski's approximation we plot in Fig the damping rate for 
both values of r] as well as the damping rate used in 0, ^] and compare with rrimag ■ One can 
observe the poor validity of the approximation for small values of g. The approximation is well 
justified if one uses the expression (p.l7|) in the range 1 < g < 2.5. 



3 Production rates 

Consider a quark-gluon plasma in which the gluons and the light quarks {u, d) are in thermal 
and chemical equilibrium. The strange quarks too are in thermal equilibrium but away from 
chemical equilibrium having very large and negative chemical potential n = Hs = Hs- The 
chemical reactions 

g + g — > s + s, (3.1) 

g + g^s + s, (3.2) 

g^s + s, (3.3) 

will then take place until chemical equilibrium is reached. The total production rate due to 
( p.lfj3.3| ), including the reversed processes, is given by @] 

SR=(l- e^^niRm^ss + Rgg^s-s + Rg^ss) (3.4) 

where 

_ [ d^Pg d^Pg d^Ps d^Ps ^o^^4x^p , p p p^ 
Kgg^s, - y (27r)32E,(27r)32E,-(27r)32E,(27r)32E,^ ^ ^ ' 

xfFD{Eg)fFD{Eg){l - /FD(i?.))(l - fFD{E,)) \M{qq ^ ss)\\ (3.5) 
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1 / d^Pl d^P2 d?Vs d^Ps n4 r/ JD TD TD TD 

^ 2 / (2vr)32Ei (27r)32E2 (27r)32E, (27r)32E,^ + P2 - P. - P. 



X 



/Bi?(i5l)/Bi?(i?2)(l - /FD(i?.))(l - /FD(i?s-)) E ^ (3-6) 



and 



~ i (27r)32E,(27r)32E,(27r)32E,^^ ' 

xfBE{Eg){l - fpDiEM^ - fFD{E,)) E \M{g ss)\^ (3.7) 

In our notation the four-momenta are denoted by capitals P, Q etc. The summation in the 
above equations extends over all colors and polarizations of the gluons and the final state 
quark-antiquark pair. We have included the Pauli blocking factors although, as long as the 
density of strange quarks is well below one, Pauli blocking does not play a significant role. 
Following Matsui, McLerran and Svetitsky we can investigate the evolution process in terms 
of the relaxation time determined near equilibrium. Therefore we have to evaluate the rates 
( P.5[j3.71 ) at /i = 0, i. e. when quarks are in both thermal and chemical equilibrium. Thermal 
field theory calculations becomes fully legitimate in this way. In particular, we can use the 
thermal quark mass given by 

ml{T) = ml{0) + ^ (3.8) 

and the thermal gluon mass ( p.l3|) . 

The processes of gluon fusion and quark-antiquark annihilation have been discussed in ref 
Q and we shall use their expressions for ( p.5| , |3.6|) . The thermal gluon decay, also discussed by 
Altherr and Seibert P, 0, can be calculated similarly. We first replace the integrations over q, 
Ps and Ps by 



d^Ps d^ps 
and change variables 



d'QdiQ^-mDeiqo) (3.9) 

d'PsS{P^-m'M)d%6{Pi-ml)e{p',) (3.10) 

Q' = Ps+Ps 

P = liPs-P^)- (3.11) 



After trivially eliminating integrals over d^Q' and d^p we find 



X 



/ dpoil - fFDiho+Pomi - fFD{\qo-Po) E Wia ^ ss)\^ (3.12) 



where the integration space is restricted by the following kinematical constraints: 



go>2m„ 0<g< (go'-4m2)i/2, ^2 ^ h _ ^ (3.I3) 
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It immediately follows that Rg^ss = if < 2ms. At the relevant temperatures the 
thermal gluon mass is not high enough to allow for decay into a strange quark pair. It is 
only because of its width that the gluon can decay. To take this into account the 5-function is 
replaced by a function, similar to the Breit-Wigner resonance. In the case of a narrow resonance 
the width of the resonance is related to the imaginary shift of the pole in the propagator in the 
complex go plane 

1 ^ 1 i2Vg^ + m^7 



which in the limit 7^0 yields the standard free particle propagator 

^ ^ , +i7rS{Q^-m^). (3.15) 



Q2 _ ^2 _ — Tn? 

Thus for a Breit-Wigner resonance with width F = 7/2 the mass-shell (5-function should be 
replaced by 

^^Q" - - - in^ %\77 -^2 ■ (3-16) 

This simple prescription cannot be directly applied to the case of thermal gluons because the 
location of the pole is determined by complicated dispersion relations ( |2.5D for transverse (T) 
and longitudinal (L) gluons. Due to instead of ( |3.16|) , we use 

r/^2_ 2^ 1 Res(AT,L)~WT,L ^ „^ 

""^^ vr (g^ _ Ren^,^)2 + Res(A^,,)-27^/ ^'-^'^ 
The matrix element is simply given by 

M{g ss) = gt,{C)u{Ps)r\av{P,), (3.18) 

where e^(C) is the polarization vector of the decaying gluon and Xa are the SU(3) matrices. 
Summing over colors and all polarizations of the gluon leads to 

Y,\M{g^ss)\^ = -VTr[(P, + m,)7^(P,-m,)7'^] 

= 16/(2m^ + g2). (3.19) 

Since the frame of the quark-gluon plasma introduces a preferred direction, it is furthermore 
necessary to distinguish between the transverse and the longitudinal components of the gluons. 
If the sum is taken over transverse or longitudinal polarization only, we find 



Y^\M{g ss)\^ = 8/[4m^ + g2(l + 44)], (3.20) 
Y^\M{g ^ srs)\' = 8g'Q\l-A%. (3.21) 
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By making use of ( glTD and (lOOl , ^ ) we find from ( ^ ) 



R 



g^ss 



27r^ 



2ms 



dqofBEiqo 



dq q 



1,11--- 



■dpo 



X(l - /FD(^go +Po))(l - fFD{^qo-Po)) 



Res(AT)-i7r 



(g2 - RenT)2 + Res(AT)-27^ 



[Ami + Q\l + A^)], 



(3.22) 



and a similar expression for Rg^ss ■ The production rate due to the gluon decay is given by 
the sum 



R 



■g^ss 



(3.23) 

If we neglect the Pauli blocking factors the integral over po can be done explicitly, leading 



g— »ss ' -^g^ss' 



to 



R: 



g^ss 



37r4 72, 



X 



dqofBEiqo) 
Res(Ar)' 



^"''(Q^ + 2m^) 



7t 



(g2 _ RenT)2 + Res(Ar)-27| 



(3.24) 



and a similar expression for 



'g^ss 



We use the full high temperature expressions for Renr,L given by ( p.l4| , prT5D and numeri- 
cally solve the dispersion relations ( |2.5| ) in order to determine Res(A2-2.) from ( ?!^ , ^.1(J| ). The 
temperature dependent gluon mass is given by (|2.13|) and the damping rate 7 = 77- 7^ is 
estimated using ( p.l7| ). 



4 Results and Conclusion 



The rates for diferent processes are depicted in Fig |^ and Fig ^. Our numerical calculation of 
the thermal gluon decay is done using equations ( |3.22| , |3.23D with ( |2.3| , |2.4| , |2^J2.1CI| ) and (|2.17| ). 
The rates for quark-antiquark annihilation and gluon fusion we calculate by making use of 
equations (3.23-3.25) in ref We fix the QCD coupling constant at the value g = 2 because 
the temperature during the time evolution is almost constant and the runing coupling effect 
is negligable. In Fig ^ the quark mass is kept fixed while in Fig |^ it varies with temperature 
according to (|3.8|). In both cases we find that the gluon fusion together with the quark- 
antiquark annihilation, dominates almost everywhere. The gluon decay process is as large as 
the gluon fusion in the narrow region around ms(0)/T = 1 only if we choose the optimistic [[7| 
parameterization of the damping rate. 

It has been shown that the time dependence of the strange-quark density can, to a great 
degree of accuracy, be described by the approximate evolution equation Q 

(4.1) 



(4.2) 



ns{t) = n^'^tanh h const 



2r 



where the relaxation time is defined as 



1 dug 

2PR 'ajJ 
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with 



The derivative of Ug with respect to /i at fixed energy density is given by 



On, 



where 



and 



= 2N 




15 



(27 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



All the quantities in ( [4.2| ) are to be evaluated at /i = 0. In Fig ^ we plot the relaxation time for 
the saturation of the strange-quark density for the massive quarks with the zero temperature 
mass ms{0) = 0.2 GeV along with the classical approximation. In this approximation Pauli 
blocking factors {l — fpo) are eliminated and the remaining Fermi-Dirac and Bose-Einstein 
distributions are replaced by the Boltzmann distribution. For comparison we also plot the 
relaxation time for the massless quarks. 

We comment here on the various approximations made in the gluon decay calculation. First 
of all, use was made of the Braaten-Pisarski resummation scheme. This is strictly valid only 
when gT « T which is clearly not the case here. This is the case with most applications 
of QCD at finite temperatures. Secondly the magnetic mass has been introduced although 
only very limited knowledge is available. It has been used to calculate the damping rate of 
a thermal gluon inside a plasma. In comparison to the calculations of references we 
keep the standard form of the Breit-Wigner distribution. The consequences of this is that the 
rate for high masses is reduced, while for low masses it is enhanced. We also avoid a rather 
heuristic assumption that the thermal quark mass is generated by gluons only. Since our rates 
are defined near equilibrium our thermal mass includes both thermal gluon and thermal quark 
contribution. Finally, a more accurate calculation of the gluon fusion process shows that the 
high mass approximation used in underestimates its contribution. 

Our main point has been that even with the parameters chosen in Ref.[^, ^ we do not 
support the claim that the gluon decay process is the dominant mechanism for strange quark 
production inside a quark-gluon plasma. To the best of our knowledge, the gluon fusion mech- 
anism is the leading process. 
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gure Captions: 

1. The damping rate and the magnetic mass versus the coupling constant. 

2. The quark production rate for thermal gluon decay for different damping rates (short 
dashed, long dashed and dot-dashed lines correspond to the damping rates depicted in 
Fig |l|) compared to the production rate for gluon fusion and quark antiquark annihilation 
(solid line). The mass rris is temperature independent. 

3. Same as Fig|^ with the thermal mass rris given by ( |3.8| ). 

4. Relaxation times for the density of massive (solid line) and massless (long dashed line) 
quarks. Corresponding relaxation times in the classical approximation are plotted with 
dashed and dotted lines respectively. 
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